Nonadiabatic electron heat pump 
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We investigate a mechanism for extracting heat from metallic conductors based on the energy- 
selective transmission of electrons through a spatially asymmetric resonant structure subject to ac 
driving. This quantum refrigerator can operate at zero net electronic current as it replaces hot by 
cold electrons through two energetically symmetric inelastic channels. We present numerical results 
for a specific heterostructure and discuss general trends. We also explore the conditions under which 
the cooling rate may approach the ultimate limit given by the quantum of cooling power. 
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I. INTRODUCTION 

The increasing miniaturization of electronic devices 
requires a deep understanding of the generation and 
flow of heat accompanying electron motion, 1,2 The quan- 
tum of thermal conductance, which is independent of 
the carrier statistics^ has been recently measured for 
phonons- and photons. 5 A practical and fundamental is- 
sue is the identification of possible cooling mechanisms 
for electron systems, a subject less developed than its 
atom counterpart.— Adiabatic electron^ and molecular 9 
pumps may provide reversible heat engines which would 
cool with minimum work expenditure. It has also been 
proposed that normal-superconductor interfaces can ef- 
ficiently cool the normal metal under appropriate condi- 
tions of electron flow. 10 ' 11 

In this paper we explore an alternative electron cool- 
ing mechanism that can operate at zero electric current 
because it relies on the idea of replacing hot electrons 
by cold electrons. The cooling concept is schematically 
depicted in Fig. [TJ An asymmetric resonant-tunneling 
structure is formed by two wells each of which hosts 
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FIG. 1: Asymmetric double- well heterostructure used for elec- 
tron ac transport calculations. Energy levels are symmetri- 
cally placed around the common Fermi level. Dominant trans- 
mission processes contributing to cooling are shown: in lead 
R hot electrons are replaced by cold electrons, all within a 
range ~ ksT around fi. 



two quasibound states. The four levels are symmetri- 
cally disposed so that the energy difference is smaller 
in the right (R) than in the left (L) well. On the 
other hand, the difference between the two upper lev- 
els is taken to be the same as that between the two 
lower ones, both being equal to the driving frequency: 
E2L - E 2 r = Em — Ei L = Ml > 0. In those condi- 
tions, electron transport is dominated by two processes: 
(i) electrons in the R electrode with energy E^r are in- 
elastically transmitted to the L electrode, where they en- 
ter with energy E 2 l = E 2 r + Ml, and (ii) electrons in 
the left with energy Em are transmitted to the right 
while also absorbing a photon. Unlike in thermionic 
refrigeration ) 12 i 13 we may assume a common chemical 
potential /1 = /il = MR- Then in the right lead one is 
effectively replacing hot electrons (with energy e > /1) by 
cold electrons (e < According to this principle, the 
right electrode is being cooled at the expense of heating 
the left electrode. This mechanism, which relies on the 
properties of coherent electron transport, may be viewed 
as the basis of a quantum refrigerator.— Under suitable 
conditions the two dominant transport mechanisms may 
cancel each other yielding a vanishing electric current, 
which prevents electrode charging. 



II. HEAT PUMP 

The classification of electrons as hot or cold depending 
on the whether its energy is above or below the chemical 
potential in its electrode is based on the property that the 
entropy variation in an infinitesimal process is given by 
TdS = dU — fidN. For independent electrons, this trans- 
lates into TdS = (e — n)dN, where e is the energy of the 
electrons being added (dN > 0) or removed (dN < 0). In 
a transport context, the entropy and temperature varia- 
tion rates are determined by the many electron scatter- 
ing processes continuously taking place at the interface. 
We always refer to the equilibrium entropy eventually 
reached in the reservoir for the new values of the con- 
served quantities energy and particle number. 

Since we are ultimately more interested in reducing 
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the temperature than the entropy, it is important to 
note that their variations are not necessarily propor- 
tional to each other. One finds CydT = (e — <j)dN, 
where Cy is the heat capacity and a = fi — T (d[i/dT) n , 
with n the particle density. In the most interesting case 
where the total electron number remains invariant on 
average (N = 0), the total entropy and temperature 
variations are proportional to each other. In the fol- 
lowing we present results for the rate of entropy varia- 
tion, knowing that it amounts to temperature variation 
in the most interesting case of constant electron number. 
Specifically, we compute the heat production rate in lead 
£ = L,RiSi!&lL!S: 



(1) 



Ni q and e q being the electron number and energy of state 
q in electrode I of chemical potential fij> = /i. 

Our goal is to understand the ac thermal transport 
properties of quantum-well heterostructures where the 
electron potential in the perpendicular z direction has 
the piecewise constant form shown in Fig. [T] while it is 
uniform in the parallel xy plane. In such a delocalized 
system, the independent-electron approximation is gen- 
erally adequate. The bottom of the right well oscillates 
as V = Vo + V ac cos (fit) while the left well operates in 
phase opposition with the same amplitude and frequency. 
To better focus on the main physical aspects, we analyze 
first transport through a single channel, later discussing 
the effect of many channels. 

Electron transport properties can be described in 
terms of scattering probabilities. Within a single-channel 
picture, the electric current flowing into lead R under ac 
driving is given byi£ 
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T$(e)h(e)-T£>(e)f R (e) , (2) 
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where fe(s) is the Fermi distribution in lead £ and T&i (e) 
is the probability for an electron to be transmitted from 
lead £' to lead £ while its energy changes from e to E+kffll, 
k being an integer number. Likewise, it can be shown 
that Eq. (TJ leads to 



k— — oo 



(3) 



+ (e + km- m)T$(e)f h (e) 
+ kMlR^(e)f R (e)] , 

where is the probability that an electron is re- 

flected in lead R from energy s to e + kMl. Invoking 
time-reversal symmetry and the monotonicity of /r(e), it 
can be proved that inelastic reflection always contributes 
to heating. Therefore any possible refrigeration of lead 
R relying on the transmission scheme depicted in Fig. [1] 
must be efficient enough to overcome the heating due 
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FIG. 2: (Color online) (a) Heat production rate in lead R for 
the structure of Fig. [T] as a function of the driving amplitude 
eFac for various lead temperatures and hfl = 1.94 meV. See 
main text for details, (b) For T = 20 K, heating contribu- 
tion from inelastically reflected electrons (solid line) and total 
electric current (line). 



to inelastic reflection. The electron scattering probabil- 
ities are calculated exactly following the transfer-matrix 
method.— 

In Fig. [2] we present numerical results for the heat pro- 
duction rate at lead R. The well lengths are 40 and 80 
nm; the heights of the barriers are Vl = Vr — 60 meV 
and Vc — 30 meV, measured with respect to the bot- 
tom of the conduction band, and their widths are 4 and 
5 nm, respectively; the difference between the bottoms 
of the two wells is Vb = 1.5 meV; the effective electron 
mass is m* — 0.07m e . This results in E2K — Em = 3.4 
meV, as determined e.g. by the dc transmission char- 
acteristics. The structure parameters have been chosen 
such that Ml = 1.94 meV coincides with E^l — ^2R and 
Em. — Em. We take // to lie half way between _Eir 
and i?2R- Clearly, the most negative heat production 
occurs for eV^c/Kl ~ 0.2. This results from a com- 
bination of nonlinearity, which yields a V a 2 c dependence 
for small F ac , and the increase of reflection heating (see 
Fig. [2t>) reinforced by the suppression of electron trans- 
mission through the dominant single-photon channels as 
eVac/Ml approaches the first zero of the first-order Bessel 
function^ The result is that \Qr\ goes through a maxi- 
mum for a moderate value of eV ac /hQ. 

Another interesting feature is that, as a function of T, 
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FIG. 3: (Color online) Cooling rate of the R electrode as a 
function of its temperature Tr for several values of Tl > Tr, 
with fi adjusted to yield zero electric current, for HQ — 1.94 
meV and eV„ c /hn = 0.2. 



the cooling rate is maximized for T ~ 20 K, which is 
roughly (E2R — Tir)/2. If the temperature is too low, 
the level 2R is empty and 1R is full, which inhibits the 
exchange of electrons. If it is too high, the cooling rate 
saturates as Tr increases, and even decreases slightly be- 
cause Tl (here equal to Tr) also increases. Later we 
argue more generally that cooling is optimized when, not 
only (E 2 r - E m )/2 but also T/2 (the half-width of the 
transmitting channels) is of order A^Tr. Here T/2 ~ 0.2 
meV, noticeably smaller than &bTr. 

The potentially most interesting scenario is that where 
cooling takes place while the net electric current is zero 
(in a classical context see Ref. [22l ). That this is not gen- 
erally the case can be inferred from the inset of Fig. [2] 
If we fix the structure and driving parameters, then the 
chemical potential and the temperatures are left as the 
independent variables. If /z is adjusted to satisfy the con- 
straint iVR = 0, the cooling rate Qr becomes a unique 
function of Tl and Tr. Figure [3] shows the resulting 
cooling rate as a function of Tr for several values of Tl. 
Remarkably, we observe that the heat production in R 
can be negative even for Tl > Tr. We conclude that it is 
technically possible to extract heat from the cold reser- 
voir and pump it to the hot reservoir with a vanishing 
net electric current. Thermodynamically, such a refriger- 
ation process requires external work, which here amounts 
to 2MI per useful scattering event and is provided by the 
classical ac source. In practice, inelastic reflection will 
further reduce the efficiency. 

In a 3D context one must generalize Eq. (3) to in- 
clude a sum over transverse modes while replacing /ir 
by A*r — h 2 K 2 /2m*, where each channel is characterized 
by its parallel wave vector k. For fixed Tl and Tr, Qr 
remains negative within a finite range of /1 values (not 
shown). This suggests that, after summing the contribu- 
tions from the many transverse channels, global cooling 
is still possible in a suitably designed 3D interface. 



III. THE QUANTUM LIMIT 

Once we have proved that it is in principle possible 
to pump heat from a cold to a hot reservoir by coher- 
ent control of electron transmission, it is natural to ask 
whether there is any fundamental limit to the maximum 
cooling rate per quantum channel which would play a role 
analogous to the quantum of electric or thermal conduc- 
tance (e 2 /h and 7r 2 fcgT/3/i, respectively). It seems evi- 
dent that the maximum cooling rate should be achieved 
in an ideal setup where a metal at temperature T is con- 
nected through a totally transparent interface to another 
metal at the same chemical potential but at zero temper- 
ature. The result is the quantum of cooling power: 

2 r°° Tr 2 k 2 T 2 

C Q = \QU, = ~j Q d sef { s) = -^-, (4) 

where /(e) = [exp(£/fc B T) + 1]" 1 and 7r 2 fc|/6/i = 473 fW 
KT 2 . Following information theory arguments, a similar 
result can be derive d 23 ' 24 . Differentiation of ([4]) yields the 
quantum of thermal conductance. Invoking only time- 
reversal symmetry and unitarity, Eq. ^ can be shown 
to satisfy (with Tr = T) 

Qr > -C Q , (5) 

for arbitrary electrodes (including fi^ ^ /xr) and driving 
parameters, thus confirming that Cq is an upper bound 
to the cooling rate. 

The quantum limit may be intuitively understood as 
follows: &)bT is the maximum amount of heat that can 
be carried away in an elementary process. Such pro- 
cesses take place at a rate ~ \Q\/ksT, which cannot ex- 
ceed h/k^T if one is to avoid effective heating caused 
by energy uncertainty. This results in \Q\ < k^T 2 /h, 
as given more precisely in ([4]). This argument suggests 
that fcfjT 2 jh is also a quantum limit for the cooling rate 
per active degree of freedom (with characteristic energy 
scale <C /crT) when cooling acts on the volume instead 
of through the surface, as e.g. in laser coolingi^^ 

The question naturally arises of whether in a setup like 
that of Fig.[T]it is possible to approach the quantum limit. 
This problem can be explored analytically within a sim- 
ple model. We neglect reflection heating and assume that 
electron transmission is dominated by two one-photon 
inelastic channels, or pipelines^ (see Fig. [IJ, named 
r u (e) = T«( £ ) and r d (e) = T^\ e ) = T^(e - Ml), 
which peak at energies T2R and Tir , respectively, always 
satisfying the unitarity requirement t u + t,j < 1. We take 
the energy origin at the middle point (Tjr + T2r)/2, 
so that T2R = — Tir = £0 > 0. Electrons entering the 
scattering region from R with initial energy ±£0 will be 
transmitted with final energy £ ± Ml through the upper 
(lower) channel. If we assume the pipelines to be sym- 
metric, r u (e) — Td(—e) = t(e), and /z = 0, we obtain 
N R = and 

Qr = -~ f dee [f R (e) - / L (e + MI)]t(s) . (6) 
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We note that at zero temperature Eq. ([6]) only yields 
heating, as should be expected. 

If we take r(e) to be a Lorentzian of width T centered 
around eo, some complications arise due to its slow decay 
for large \e — £q\. For instance, for large enough f2, we 
always find heating Qr oc lnf2. On the other hand, for 
small ft, Qr < if and only if Tl < Tr. We conclude 
that, in the interesting case Tl > Tr, cooling of the R 
electrode can only occur within a finite range of £1 values. 
This range shrinks to zero for Tl large. 

An interesting question is whether, given two elec- 
trodes with Tl > Tr, it is always possible to design an 
ac resonance structure yielding Qr < 0, and whether 
|Qr| can ever approach the quantum limit. In Eq. (|6j). 
g{e) = e [fn(e) — + ti&l)] > only in the interval 
< e < e = ftQT R /(T L - T R ). For T L -> T R , we have 
e — > oo; however, the integrand decays exponentially on a 
scale ~ kftT after having peaked at e ~ fceT. Therefore, 
cooling comes effectively from the interval < e < £\, 
where ex = min{e, 2fc B T} and f = (T L + T R )/2. To 
potentiate the contribution from that segment, we may 
design r(e) to be centered at £o — £i/2. If V — *■ 0, Qr 
is guaranteed to become negative, although with a van- 
ishing magnitude |Qr| oc T. A typical optimal value is 
r ~ e\- We conclude that the cooling rate is maximized 
for £ ~ T/2 - ex/2. The peak at k B T ~ e for T L = T R 
is confirmed by the lowest curve of Fig. [3) 

If t(s) decays sufficiently fast away from the re- 
gion where g(e) > 0, one may estimate \Qn\ ~ 
-(2//i)eiT max 5 max . For Tl — > T R we have both ex and 
ffmax of order fceT, assuming Ml 3> k^T . In those con- 
ditions, |Qr| ~ Cq provided r max is close to unity. By 
contrast, the cooling rate cannot approach the quantum 
limit if Tl grows substantially above Tr or if r(e) decays 
slowly, like in a Lorentzian resonance, since then the con- 
tribution from g(e) < cannot be neglected. 

Careful inspection of Eq. ^ reveals that Qr increases 
monotonically as Tr decreases. Thus, if we start cooling 
the R electrode, Qr begins to increase until it eventually 
becomes zero. At that point, no further cooling is possi- 
ble. We have reached the lowest possible temperature for 
the refrigeration process defined by r(e) and ft. We are 
limited by the lack of sufficient energy resolution: when 
/cbTr becomes small compared with the linewidth T, no 
heat pumping is possible for Tl > Tr. 

IV. DISCUSSION 

The quantum refrigerator which we have investigated 
may be viewed as a realization of Maxwell's demor*2& as 



it selectively lets hot electrons out while it only lets cold 
electrons in. The required work is provided by the exter- 
nal ac source which, combined with the spatial asymme- 
try of the structure, rectifies electron motion. The work 
might also be extracted from a hot ohmic resistor^ Al- 
ternative schemes to that of Fig. Q] are of course possible: 
One may design two superlattices each of them having 
two narrow bandwidths yielding a similar level distri- 
bution. A potential advantage of such a device would 
be that, away from resonance, transmission would de- 
cay fast. Thus it would show interesting features such as 
cooling for arbitrary large f2 and, as discussed above, the 
guaranteed existence of a driving structure that brings 
the cooling rate close to the quantum limit. A non- 
resonant mechanical mismatch at the interface would pre- 
vent phonons from short-circuiting electron cooling dur- 
ing operation close to such a limit. 

In conclusion, we have identified a mechanism for nona- 
diabatically pumping heat from a cold to a hot electron 
reservoir which is based on the coherent control of elec- 
tron ac transport and which can operate at zero aver- 
age electric current. On the basis of electron transport 
considerations, the quantum of cooling power Cq has 
been shown to be an upper bound to the cooling rate 
per quantum channel. We have investigated the case of 
Lorentzian resonances, where approaching the quantum 
limit is generally not feasible. We have noted however 
that with sharper resonances it is always possible to de- 
sign a driven interface that provides cooling at a rate 
close to the quantum limit. 
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